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Efforts  to  harness  wind  energy  on  a  large  scale  have  gained  momentum  across  the  world.  By  the  end  of 
December  2013,  a  cumulative  capacity  of  more  than  300  GW  of  wind  power  projects  had  been  installed 
all  over  the  world.  One  of  the  key  aspects  involved  in  implementing  wind  power  projects  is  the  analysis 
of  wind  speeds  distributions  observed  or  recorded  and  assessment  of  annual  energy  output  from  the 
wind  turbines.  The  wind  speed  frequency  distribution  is  generally  assumed  to  follow  two-parameter 
Weibull  Distribution,  in  general,  across  the  world,  annual  energy  generation  estimations  of  a  wind 
turbine  at  a  given  site  are  assessed  on  the  basis  of  Weibull  Distribution.  However,  in  this  paper,  based  on 
a  robust  analysis  carried  out  on  over  208  measurement  sites  in  India,  we  show  that  multi-peak  Gaussian 
distribution  functions  are  a  significantly  improved  representation  of  observed  wind  speed  distributions. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  worldwide  installed  capacity  of  wind  power  was  a  little 
more  than  300  GW  at  the  end  of  December  2013  [1],  The 
electricity  generated  from  these  windfarms  accounts  for  more 
than  3.5%  of  the  global  electricity  consumption  [1  ].  Geographically, 
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there  is  now  a  widespread  deployment  of  modern  wind  turbine 
technology  across  the  globe  and  nearly  all  the  countries  are 
including  wind  energy  in  their  plans  and  policies.  Over  the  last 
5-10  years,  utility-scale  wind  electricity  generation  has  emerged 
as  a  mature  mainstream  energy  technology. 

A  prerequisite  for  setting  up  a  windfarm  is  detailed  wind 
resource  assessment  at  the  site,  which  includes  measurement  of 
wind  speeds  and  other  climatic  conditions  for  a  minimum  3-year 
period.  However,  given  the  pace  at  which  windfarms  are  being 
planned  and  the  commensurate  speed  with  which  policies,  pro¬ 
grams  and  the  legal,  statutory  and  regulatory  frameworks  should 
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be  set  in  place,  the  stakeholders  and  agencies  involved  are  always 
in  a  hurry  and  often  try  to  reduce  the  development  time.  An 
important  component  in  the  development  time  is  the  wind  speed 
measurement  period.  Reduction  in  measurement  period  can  be 
achieved  by  means  of  mathematical  modelling  of  different  aspects 
of  wind  resource.  Some  of  these  aspects  are  vertical  and  horizontal 
extrapolation  of  wind  regimes,  regional  assessment  of  large 
geographical  areas,  wind  distribution  modelling,  correlations  with 
nearby  measurements,  etc.  With  these  techniques,  it  is  possible  to 
reduce  not  only  the  measurement  but  also  unnecessary  expenses. 

Here,  mathematical  modelling  of  wind  speed  distributions 
plays  an  important  role.  Assessments  of  annual  energy  output 
(AEO)  from  wind  turbines  are  often  made  using  mathematical 
models  in  which  a  mathematical  function  is  assumed  to  represent 
the  actual  wind  speed  frequency  distribution.  This  enables  AEO 
computations  to  be  made  by  using  only  a  few  parameters  and  the 
mathematical  function  mimics  the  frequency  distribution  of  wind 
speeds. 

Though  there  are  different  kinds  of  functions  that  could  be 
fitted  to  wind  speeds  measured  at  a  site,  the  two-parameter 
Weibull  distribution  has  wide  acceptance  across  the  world.  In 
spite  of  scientific  improvements,  standardization  of  measurement 
procedures  such  as  International  Electrotechnical  Commission 
(1EC)  [2],  variance  in  AEO  assessments  from  actual  output  of 
windfarms  continues  to  be  an  area  of  major  concern  in  the  wind 
industry.  The  purpose  of  the  work  reported  here  is  to  assess  the 
suitability  of  multi-peak  Gaussian  function  to  wind  speed  distri¬ 
butions  in  place  of  Weibull  distribution  function  and  to  compare 
the  two. 

Our  analysis  with  wind  speed  data  from  208  locations,  spread 
all  over  India,  indicates  that  wind  speed  distributions  have  a  much 
better  fit  with  Gaussian  multi-peak  functions  as  compared  to 
Weibull  distribution. 


2.  A  review  of  different  models 

Justus  et  al.  [3]  have  justified  the  use  of  Weibull  distribution 
and  its  subset,  the  Rayleigh  distribution,  which  is  the  Weibull 
distribution  with  “k= 2".  Innumerable  books  and  publications 
[4-7]  have  used,  projected  or  justified  Weibull  distribution  or 
Rayleigh  distributions  as  good  representation  of  wind  speed 
distributions.  Sedefian  [8]  has  assumed  Weibull  distribution  in 
his  presentation  of  methods  of  extrapolating  wind  distributions  to 
different  heights.  Pang  et  al.  [9]  have  reiterated  the  reasons  for 
using  Weibull  distribution  given  by  Justus  et  al.  [3],  Pang  et  al. 
have  estimated  three  parameters  of  Weibull  distribution  using 
Markov  Chain  Monte  Carlo  (MCMC)  method  and  Max  Likelihood 
(ML)  methods.  They  have  concluded  that  for  the  data  analysed  by 
them,  there  was  little  evidence  that  the  three-parameter  Weibull 
model  was  necessary,  thus  implying  that  the  two-parameter 
Weibull  distribution  described  above  was  adequate  for  approx¬ 
imating  wind  speed  distribution. 

Both  Lysen  and  Justus  have  also  mentioned  Rayleigh  distribu¬ 
tion  which  is  a  special  case  of  Weibull  distribution  with  k= 2.  The 
Rayleigh  distribution  offers  additional  advantage  over  Weibull 
distribution  due  to  the  fact  that  the  computations  become  more 
simplified  as  only  one  variable,  i.e.,  ‘mean  wind  speed’,  has  to  be 
used  in  the  computations.  This  could  have  a  significant  implication 
on  computational  time  and  effort  when  assessments  are  being 
made  for  a  geographical  area  or  at  a  regional  level. 

A  study  on  Antarctica  has  used  Rayleigh  distribution  [10],  In  the 
assessment  of  California  offshore  assessment  [11]  and  the  Global 
wind  potential  assessment  also,  Rayleigh  distribution  has  been 
assumed. 


In  more  recent  work  by  Celik  and  Muneer  [12],  a  critical 
evaluation  of  different  frequency  distributions  has  been  carried 
out.  In  addition  to  two-parameter  Weibull  and  Rayleigh  distribu¬ 
tions  they  have  also  evaluated  the  applicability  of  three-parameter 
Weibull  distribution,  lognormal  distribution  and  bimodal  Weibull 
distribution.  For  the  data  analyzed,  they  have  found  bimodal 
Weibull  distribution  provides  the  best  fit.  However,  the  short¬ 
coming  with  this  work  is  that  it  uses  data  from  only  one  site  and 
secondly  the  objective  of  the  paper  seems  to  be  development  of  a 
method  of  arriving  at  a  score  to  figure  out  which  model  is  best 
suited  to  a  given  wind  distribution  and  not  to  establish  as  to  which 
model,  in  general  and  for  wider  applicability,  is  the  appropriate 
approximation  for  representing  wind  speed  distributions. 
Celik  [13]  assessed  the  error  in  energy  estimation  of  small  wind 
power  systems  using  Weibull  distribution  to  be  of  the  order  of 
2.79%.  Kollu  et  al.  [14]  have  evaluated  three  mixture  probability 
density  functions  Weibull-extreme  value  distribution  (GEV),  Wei- 
bull-lognormal,  and  GEV-lognormal  distributions  for  their  suit¬ 
ability  to  wind  distributions  and  concluded  that  mixture 
distributions  are  able  to  provide  better  fit.  Akdag  et  al.  [15]  have 
evaluated  two-parameter  Weibull  distribution  and  two- 
component  mixture  Weibull  distribution  with  five  parameters 
and  concluded  a  mixture  of  two  Weibull  distributions  is  more 
suitable  for  the  description  of  such  wind  conditions  and  could 
offer  less  relative  errors  in  determining  the  annual  mean  wind 
power  density. 

It  is  interesting  to  note  that  a  mixture  of  distributions  turns  out 
to  be  a  better  fit  than  a  single  Weibull  distribution.  It  is  a  known 
fact  that  wind  speed  distributions  are  often  multi-modal  or  bi¬ 
modal  particularly.  Such  bi-modal  or  multipeak  distributions  will 
always  present  difficulties  in  fitting  them  to  a  single  Weibull 
distribution. 

With  the  exception  of  more  recent  work  of  Celik  and  Muneer 
[12],  Kollu  et  al.  [14]  and  Akdag  et  al.  [15],  who  have  either  carried 
out  a  critical  evaluation  of  different  frequency  distributions  w.r.t. 
wind  energy  or  have  attempted  mixtures  of  two  or  more  prob¬ 
ability  distribution  functions,  literature  search  carried  out  so  far 
continues  to  point  overwhelmingly  towards  two-parameter 
Weibull  distribution  or  the  Rayleigh  distribution  as  the  widely 
accepted  choice  for  representing  wind  speed  distribution.  How¬ 
ever,  in  our  analysis  presented  here,  we  find  that  multipeak 
Gaussian  distributions  are  a  significantly  better  fit  than  Weibull 
and  Rayleigh  distributions. 


3.  Theory  and  formulations 

3.1.  Probability  density  function 

The  probability  density  function  (PDF)  is  a  function  that  can  be 
integrated  to  obtain  the  probability  of  occurrence  of  the  variable 
having  a  value  in  a  given  class  interval.  The  probability  density 
function  is  nonnegative  everywhere,  and  its  integral  over  the 
entire  space  is  equal  to  one. 

A  probability  density  function  is  most  commonly  associated 
with  absolutely  continuous  univariate  distributions.  A  random 
variable  X  has  density  f  where  /  is  a  non-negative  function 
[16-18],  if 

rb 

P[a<X<b]=  J  f(x)dx  (1) 

Hence,  if  F  is  the  cumulative  distribution  function  of  “x”,  then 
F(x)=  [  f(u)du  (2) 
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and  (if/ is  continuous  at  “x”) 

/(*)  =  ^F(X)  (3) 

Intuitively,  one  can  think  of  f[x )  dx  as  being  the  probability  of 
“x”  falling  within  the  infinitesimal  interval  [x,  x+dx[. 

3.2.  Normal  distribution 

A  random  variable  which  has  a  normal  distribution  with  a  mean 
m=0  and  a  standard  deviation  <r=  1  is  referred  to  as  Standard  Normal 
Distribution. 


3.3.  Weibull  distribution 


The  Weibull  distribution  is  widely  used  in  engineering  and 
scientific  analysis.  It  is  used  in  survival,  reliability  and  failure 
analysis,  manufacturing  and  delivery  times,  extreme  value  theory 
and  weather  forecasting.  In  wind  energy,  it  is  a  well-established 
fact  that  a  wind  speed  frequency  distribution  can  be  approximated 
by  a  Weibull  distribution  [19],  which  is  described  by  shape 
parameter  “k"  and  scale  parameter  “c”.  Rayleigh  distribution,  a 
special  case  of  Weibull  distribution  with  k=2,  is  also  widely  used 
to  represent  wind  distributions. 

The  Weibull  probability  density  function  (PDF)  for  wind  speed 
“v”  is  given  by 


f(v)  =  (k/c)(^ 


exp 


(4) 


In  the  above  (k  >  0,  v  >  0,  c  >  1 ). 

The  relationship  between  scale  parameter  “c",  shape  parameter 
“k"  and  mean  wind  speed  “um”  is  given  by 


vm  =  cr (1  +  k  ') 


(5) 


where  “T”  is  the  Gamma  function. 

The  energy  output  from  the  wind  turbine  over  time  period  “T” 
is  given  by 

rV,  rvou, 

E  =  T  P(V)f(v)dv+  TPr  f(v)dv  (6) 

Jv m  Jv, 

where  “P(V)”  is  the  power  curve  function  of  the  wind  turbine,  “Pr" 
is  its  rated  power,  “ vin ”  is  the  cut-in  wind  speed,  “v"  is  the  rated 
wind  speed  and  “vout”  is  the  cut-out  wind  speed  and  is  given 
by  Eq.  (4).  A  form  of  Eq.  (6)  shown  in  Eq.  (7)  can  be  used  to 
compute  annual  energy  output  (AEO)  from  a  given  wind  turbine: 

AEO  =  8760  £  Pi  Vi  (7) 


In  Eq.  (7),  “V"  is  the  frequency  for  different  wind  speed  class 
intervals  of  a  wind  speed  frequency  distribution  that  can  be 
derived  from  Eq.  (4)  for  frequency  distributions  with  different 
“k”  and  “c”  and  “Vm”,  n  is  the  number  of  class  intervals  and  “P,”  is 
the  value  of  power  generated  corresponding  to  the  mid  value  of 
the  ith  class  interval  and  is  obtained  from  the  wind  turbine  power 
curve,  P(V). 


3.4.  Gaussian  distribution 


Apart  from  fitting  Weibull  and  Rayleigh  distributions  to  the 
percent  frequency  distributions,  we  also  fitted  Gaussian  distribu¬ 
tion.  A  Gaussian  distribution  is  given  by 

g(x)=  ae~(x~b)2/2c2  (8) 

where  a,  b,  and  c  are  the  curve  parameters  that  control  the  shape 
of  the  function,  ‘a’  is  the  height  of  the  curve  at  the  maximum,  ‘b’  is 
the  value  of  ‘x’  where  the  maximum  occurs,  and  ‘c’  controls  the 


width  of  the  curve,  that  is,  how  rapidly  it  falls  to  zero.  To  use  the 
Gaussian  as  a  probability  distribution  function,  we  require  that  the 
total  area  under  the  curve  equals  one.  The  normalized  Gaussian 
distribution  function  is  given  by 

G(x)=  _^e-(*-rt2/2<P  (9) 

Os/2n 

The  curve  fitting  software  ORIGIN  9  (www.originlab.com)  enables 
such  a  curve  fitting  over  the  graph  plotted  on  its  graphing  tool. 

Goshtasby  and  O’Neill  [21]  have  shown  that  a  sequence  of 
uniformly  spaced  points  {x,:  i=l,...,N)  with  associated  data  values 
fyfc:  k=\,...,N)  can  be  approximated  by  a  sum  of  Gaussians  in  the 
form 

y=f(x)=  £  (10) 

i  =  1 

After  fitting  the  wind  speed  frequency  distribution  to  Gaussian 
functions  as  in  Eq.  (10),  AEO  can  be  computed  by  Eq.  (7)  as  was 
done  in  the  case  of  Weibull  distribution. 


3.5.  R-squared  value 


The  R-squared  value,  a  measure  of  goodness  of  fit  used  to 
compare  the  results  of  fitting  different  types  of  distributions  to  the 
actual  distribution,  is  given  by 


R  =  1  — 


SSE 

SST 


(11) 


where  SSE,  also  known  as  summed  square  of  residuals,  is  given  by 


SSE  =  £(  Xi-y,)2 
i 


(12) 


where  “Xj”  are  the  actual  or  observed  frequency  distribution  values 
for  N  class  intervals  and  “y"  is  the  predicted  value  (in  this  case 
through  the  Gaussian  fit). 

SST  is  the  total  sum  of  squares  given  by 

SST=  £(xi-P)2  (13) 

i 


4.  Data  and  analysis 

In  order  to  compare  suitability  of  multi-peak  Gaussian  distri¬ 
bution  with  Weibull  and  Rayleigh  distribution,  we  carried  out 
analysis  with  data  measured  at  208  sites  in  India.  The  locations  of 
these  measurement  sites  are  shown  in  the  map  in  Fig.  1.  The  data 
has  been  measured,  analysed  and  published  by  CWET.  The  pub¬ 
lications  by  CWET  present  percentage  frequency  distribution, 
shape  parameter  k,  scale  parameter  c,  mean  annual  wind  speeds 
at  different  heights  and  air  density. 

These  measurements  are  being  carried  out  under  a  National 
Program  of  the  Ministry  of  New  and  Renewable  Energy  (MNRE) 
that  involves  setting  up  wind  masts  at  varying  heights.  Till  date 
nearly  2080  [22]  wind  measurement  masts  at  20  m,  25  m  and 
50  m  have  been  set  up  and  the  data  is  being  periodically  published 
by  the  Chennai-based  Centre  for  Wind  Energy  Technologies 
(CWET),  which  is  the  scientific  wing  of  MNRE.  The  first  compila¬ 
tion  [23]  was  published  in  1983.  Over  the  years,  such  publications 
[24-31]  have  continued  to  provide  updated  processed  information 
on  wind  speeds  in  India. 

Under  this  programme,  all  the  wind  monitoring  stations  are 
equipped  with  identical  instrumentation.  Each  installation  con¬ 
sists  of  a  20/25  m  tall,  75  mm  diameter,  guyed  tubular  steel  mast 
with  booms  fixed  at  two  levels  10  m  and  20/25  m  above  ground 
and  a  wind  resource  data  logging  system.  Now  wind  masts  with 
height  up  to  120  m  above  the  ground  are  also  being  set  up.  Data  is 
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Fig.  1.  Overview  of  analysed  sites  (208)  location  on  Indian  map  published  by  Centre  of  Wind  Energy  Technology  (CWET),  Chennai. 


Fig.  2.  Comparisons  of  R-square  value  of  65  sites  based  on  Weibull,  Rayleigh  and 
Gaussian  distribution  fit. 


-•--Wb 
— 3K —  Ra 


Sites 


sampled  at  a  frequency  of  0.5  Hz  (2  s)  and  averaged  over  one  hour. 
The  data  loggers  are  set  to  store  the  wind  speed  in  ‘m/s’  (meter- 
per-second),  wind  direction  and  standard  deviation  of  wind  speed 
at  two  levels.  In  order  to  ensure  that  data  collected  is  both  reliable 
and  accurate,  the  installation,  maintenance  and  periodic  inspec¬ 
tion  of  the  instruments  and  the  collection  of  chips  have  been 
carried  out  by  CWET  under  rigidly  controlled  conditions. 

We  have  used  percentage  frequency  distribution,  shape  para¬ 
meter  k,  scale  parameter  c  and  mean  annual  wind  speeds  from  the 
above  publications  [17,  21-28]  for  the  208  sites  spread  all  over  the 
country.  The  percentage  frequency  distribution  presented  in  these 


Fig.  3.  Error  in  AEO  estimations  based  on  actual  in  comparison  with  Gauss,  Weibull 
and  Rayleigh  distribution  functions. 

publications  is  treated  as  the  actual  wind  speed  distribution  at  the 
site  and  then  the  R-squared  (R2)  values,  that  are  a  measure  of  the 
goodness  of  fit,  have  been  computed  for  fitting  three  different 
distribution  functions  to  actual  wind  speed  distribution,  i.e., 

•  Weibull 

•  Rayleigh 

•  Multi-peak  Gaussian 
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We  have  also  computed  annual  energy  output  (AEO)  from  a 
given  wind  turbine  (Suzlon  S-88)  for  these  sites  to  study  the 
sensitivity  of  AEO  with  respect  to  the  frequency  distributions 
assessed  and  also  to  assess  the  error  in  estimation  in  comparison 
with  the  assessment  made  with  actual  frequency  distribution.  The 
^-squared  values  between  actual  measurement  for  these  sites 
(percent  frequency  distribution)  and  the  multi-peak  Gauss,  Ray¬ 
leigh  and  Weibull  distributions  are  presented  in  Fig.  2  (for  better 
visibility  65  sites  are  presented).  The  Mean  Absolute  Error  in  AEO 
computations  when  compared  with  AEO  computations  from 
actual  distribution  were  found  to  be  4%  with  Gauss,  6%  with 
Weibull  and  8%  with  Rayleigh  distribution.  The  plot  of  error  in  AEO 


estimation  for  a  dataset  comprising  47  sites  in  Fig.  3  shows  that 
Gaussian  multi-peak  distribution  results  in  minimum  error  in 
energy  estimations. 


5.  Results  and  discussion 

It  can  be  seen  from  Eq.  (7)  that  AEO  computations  are  sensitive 
to  the  kind  of  distribution  assumed  for  wind  speeds. 

Fig.  4  presents  percentage  frequency  wind  speed  data  in 
decimals  at  one  of  the  sites  out  of  the  208  sites  in  Karnataka  State 
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— Ga 
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■  x  Ac 

— *  -Ra 


Wind  Speed  Class  Interval  (m/s) 


Wind  speed  class  interval  (m/s) 


Fig.  4.  General  representation  of  percent  frequency  wind  speed  distribution 
function. 


Fig.  5.  A  comparison  of  different  frequency  distribution  functions  (site: 
Kapattagudda). 


Table  1 

R2  value  computation  using  Gaussian  distribution  function  at  Channavadayanpura  using  multi-peak  technique  for  curve  fitting. 


Site:  Channavadayanpura 

Mean  wind  speed:  6.3 

Weibull  shape  factor  k:  2.4 

Guass-Actual  fit 

Output  from  ORIGIN  on  Gaussian  fit  (Fig.  5) 

SSE  (Actual-Guass) 

0.000463 

Area  (A) 

Mean  (xc) 

Width  (w) 

SSE  (Actual-Guass) 

0.058090 

Peak  1 

0.31349 

3.10299 

3.12733 

SSE/SST 

0.007964 

Peak  2 

3.25 

6.91794 

5.09 

R-squared 

0.992036 

Peak  3 

0.3704 

6.54123 

3.16709 

Wind  speed  (m/s) 

Actual  %  Freq  distribution 

Guass  distributions 

Guass  peaks  summed 

SSE  calc. 

SST  calc. 

Actual 

GS  PK  1 

GS  PK  2 

GS  PK  3 

Sum  PKs 

0.00 

0 

0.011165 

0.001266 

0.000018 

0.012450 

0.000155 

0.002776 

1.00 

0.05 

0.032376 

0.003410 

0.000205 

0.035991 

0.000196 

0.000007 

2.00 

0.063 

0.062365 

0.007874 

0.001528 

0.071767 

0.000077 

0.000106 

3.00 

0.107 

0.079808 

0.015576 

0.007656 

0.013041 

0.000016 

0.002950 

4.00 

0.118 

0.067847 

0.026405 

0.025747 

0.120000 

0.000004 

0.004266 

5.00 

0.137 

0.038317 

0.038358 

0.058110 

0.134785 

0.000005 

0.007109 

6.00 

0.148 

0.014376 

0.047747 

0.088021 

0.150144 

0.000005 

0.009085 

7.00 

0.145 

0.003583 

0.050931 

0.089480 

0.143994 

0.000001 

0.008522 

8.00 

0.109 

0.000593 

0.046555 

0.061049 

0.108194 

0.000001 

0.003171 

9.00 

0.063 

0.000065 

0.036462 

0.027953 

0.064481 

0.000002 

0.001060 

10.00 

0.033 

0.000005 

0.024472 

0.008590 

0.033067 

0.000000 

0.003870 

11.00 

0.017 

0.000000 

0.014075 

0.001772 

0.015847 

0.000001 

0.001273 

12.00 

0.007 

0.000000 

0.006936 

0.000245 

0.007182 

0.000000 

0.002087 

13.00 

0.003 

0.000000 

0.002929 

0.000023 

0.002952 

0.000000 

0.002469 

14.00 

0.001 

0.000000 

0.001060 

0.000001 

0.001061 

0.000000 

0.002671 

15.00 

0 

0.000000 

0.003290 

0.000000 

0.000329 

0.000000 

0.002776 

16.00 

0 

0.000000 

0.000087 

0.000000 

0.000087 

0.000000 

0.002776 

17.00 

0 

0.000000 

0.000020 

0.000000 

0.000020 

0.000000 

0.002776 

18.00 

0 

0.000000 

0.000004 

0.000000 

0.000004 

0.000000 

0.002776 

Mean 

0.052684 

488 


J.  Hossain  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  34  (2014)  483-490 


Table  2 

R2  value  computation  using  Weibull  and  Rayleigh  distribution  function  at  Channavadayanpura  on  multi-peak  wind  speed  distribution. 


Site:  Channavadayanpura 

Mean  wind  speed:  6.3  m/s 

Wind  speed  (m/s)  Actual  %  freq  distribution 

SST  (Actual-WEIBULL) 
SST  (Actual-WEIBULL) 
SSE/SST 
//-.squared 

Weibull 

0.0032 

0.0581 

0.0550 

0.9450 

SSE  calc. 

SST  calc. 

SST  (Actual-Rayleigh) 
SST  (Actual-Rayleigh) 
SSE/SST 

R-squared 

Rayleigh 

0.005415 

0.058090 

0.093225 

0.906775 

SSE  calc. 

SST  calc. 

0.00 

0 

0.0000 

0.0000 

0.0028 

0.0000 

0.0000 

0.0028 

1.00 

0.05 

0.0215 

0.0008 

0.0000 

0.0388 

0.0001 

0.0000 

2.00 

0.063 

0.0546 

0.0001 

0.0001 

0.0732 

0.0001 

0.0001 

3.00 

0.107 

0.0891 

0.0003 

0.0030 

0.0995 

0.0001 

0.0030 

4.00 

0.118 

0.1175 

0.0000 

0.0043 

0.1154 

0.0000 

0.0043 

5.00 

0.137 

0.1343 

0.0000 

0.0071 

0.1207 

0.0003 

0.0071 

6.00 

0.148 

0.1369 

0.0001 

0.0091 

0.1165 

0.0010 

0.0091 

7.00 

0.145 

0.1261 

0.0004 

0.0085 

0.1051 

0.0016 

0.0085 

8.00 

0.109 

0.1055 

0.0000 

0.0032 

0.0892 

0.0004 

0.0032 

9.00 

0.063 

0.0806 

0.0003 

0.0001 

0.0717 

0.0001 

0.0001 

10.00 

0.033 

0.0562 

0.0005 

0.0004 

0.0546 

0.0005 

0.0040 

11.00 

0.017 

0.0358 

0.0004 

0.0013 

0.0397 

0.0004 

0.0013 

12.00 

0.007 

0.0209 

0.0002 

0.0021 

0.0274 

0.0002 

0.0021 

13.00 

0.003 

0.0111 

0.0001 

0.0025 

0.0181 

0.0001 

0.0025 

14.00 

0.001 

0.0054 

0.0001 

0.0028 

0.0114 

0.0000 

0.0027 

15.00 

0 

0.0024 

0.0000 

0.0028 

0.0069 

0.0000 

0.0028 

16.00 

0 

0.0009 

0.0000 

0.0028 

0.0040 

0.0000 

0.0028 

17.00 

0 

0.0003 

0.0000 

0.0022 

0.0022 

0.0000 

0.0028 

18.00 

0 

0.0001 

0.0000 

0.0012 

0.0012 

0.0000 

0.0028 

Fig.  6.  Modelled  (Weibull,  Rayleigh  and  Gaussian  distributions)  fit  for  observed 
data  values  of  percentage  frequency  distribution  of  wind  speed  (m/s)  averaged 
hourly  at  Amberi  (Maharashtra). 


-  Channavadayanpura  in  coastal  Karnataka.  We  have  used  percent 
frequency  data  (which  represents  actual  distribution)  and  Weibull 
shape  parameter  values  for  208  sites  in  our  analysis.  Table  1 
presents  the  R2  value  which  has  been  computed  at  Channava¬ 
dayanpura  using  multi-peak  Gaussian  distribution  function  tech¬ 
nique  for  curve  fitting. 

A  comparison  of  the  three  different  distributions  with  actual 
percentage  frequency  distribution  at  a  site  Kapattaguda  is  shown 
in  Fig.  5. 

A  vast  study  of  percentage  frequency  distribution  of  wind 
speed  data  has  been  carried  out  and  different  types  of  peaking 
in  the  distributions  have  been  observed,  e.g.  uni-mode,  bi-mode 
and  multi-mode.  We  have  tried  to  fit  these  distributions  with 
Weibull,  Rayleigh  and  Gaussian  distribution  functions  and 


Class  Interval  (m/s) 

Fig.  7.  Modelled  (Weibull,  Rayleigh  and  Gaussian  distributions)  fit  for  observed 
data  values  of  percentage  frequency  distribution  of  wind  speed  (m/s)  averaged 
hourly  at  Ulindakonda  (Andhra  Pradesh). 


compared  the  best  fit  results.  Table  2  presents  the  R 2  values  that 
have  been  computed  at  Channavadayanpura  using  Weibull  and 
Rayleigh  distribution  functions.  In  this  research  paper  R2  values 
are  used  to  compare  the  best  fit  model  with  actual  wind 
speed  data. 

In  order  to  show  our  results  in  pictorial  view  two  sites  have 
been  selected  for  each  peak  mode. 

5.3.  Uni-mode  curve 

At  Amberi  (Maharashtra)  and  Ulindakonda  (Andhra  Pradesh) 
wind  speed  distribution  is  uni-modal  as  shown  in  Figs.  6  and  7, 
respectively.  In  both  Figs.  6  and  7,  Weibull,  Rayleigh  and  Gauss 
distributions  show  approximately  similar  R2  values  varying  from 
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0.95  to  0.99.  The  R2  value  is  closer  to  “1.0”  in  Gaussian  fit 
probability  distribution.  The  high  value  of  R2  shows  good  relation¬ 
ship  between  observed  and  modelled  values. 


Fig.  8.  Modelled  (Weibull,  Rayleigh  and  Gaussian  distributions)  fit  for  observed 
data  values  of  percentage  frequency  distribution  of  wind  speed  (m/s)  averaged 
hourly  at  Pushpathur  (Tamil  Nadu). 


Fig.  9.  Modelled  (Weibull.  Rayleigh  and  Gaussian  distributions)  fit  for  observed 
data  values  of  percentage  frequency  distribution  of  wind  speed  (m/s)  averaged 
hourly  at  Mylampati  (Tamil  Nadu). 


Fig.  10.  Modelled  (Weibull,  Rayleigh  and  Gaussian  distributions)  fit  for  observed 
data  values  of  percentage  frequency  distribution  of  wind  speed  (m/s)  averaged 
hourly  at  Aundhewadi  (Maharashtra). 


Fig.  11.  Modelled  (Weibull,  Rayleigh  and  Gaussian  distributions)  fit  for  observed 
data  values  of  percentage  frequency  distribution  of  wind  speed  (m/s)  averaged 
hourly  at  Warshamadi  (Gujarat). 


5.2.  Bi-mode  and  multi-mode  curve 

Pushpathur  (Tamil  Nadu)  and  Mylampati  (Tamil  Nadu)  sites 
have  been  selected  for  bi-modal  (two  prominent  peaks)  illustra¬ 
tion  in  Figs.  8  and  9  respectively.  Here  we  see  that  while  the  R2 
value  is  close  to  1  in  case  of  Gauss  function,  it  is  much  lesser  in 
case  of  Weibull  and  Rayleigh  distribution  functions.  We  find  a 
similar  result,  in  multi-modal  (more  than  two  prominent  peaks) 
illustration  in  Figs.  10  and  11  for  Aundhewadi  (Maharashtra)  and 
Warshamedi  (Gujarat)  sites,  respectively.  This  exercise  has  been 
carried  out  over  208  measurement  sites. 

Bi-modal  and  multi-modal  curve  fitting  with  Weibull,  Rayleigh 
and  Gauss  frequency  distribution  functions  has  been  compared. 
The  R2  value  is  found  to  be  significantly  higher  for  Gaussian  fit 
when  compared  with  Weibull  and  Rayleigh  distribution  functions. 

In  wind  resource  assessment,  determination  of  accurate  annual 
energy  output  (AEO)  is  of  great  importance.  Since  the  AEO 
determination  is  linked  to  distribution  curve  fit,  it  requires  careful 
treatment  in  analysis  and  estimation.  We  feel  that  the  most 
important  finding  from  our  work  reported  here  is  the  greater 
suitability  of  Gaussian  Multi  Peak  functions  as  compared  to 
Weibull  and  Rayleigh  functions  to  a  wide  range  of  wind  speed 
datasets,  covering  almost  the  whole  of  India.  Given  the  size  of  our 
dataset,  we  can  claim  that  this  finding  is  indeed  universally 
applicable. 


6.  Conclusions 

Our  main  finding  is  that  the  wind  distribution  function  is  better 
approximated  by  multi-peak  Gaussian  curves  in  comparison  with 
customarily  assumed  Weibull  and  Rayleigh  distribution  functions. 
This  finding  is  based  on  analysis  carried  out  with  measurements  at 
208  sites  with  varying  wind  regimes  throughout  the  Indian 
subcontinent.  The  approach  of  AEO  assessment  as  presented  in 
the  study  using  multi-peak  Gaussian  function  can  be  efficiently 
and  effectively  adopted  for  greater  accuracy.  The  robustness  of 
datasets  and  the  findings  reveal  that  the  phenomenon  is  univer¬ 
sally  applicable. 

There  is  worldwide  growth  and  expansion  in  the  deployment 
of  wind  energy  technologies  and  in  this  process  accuracy  in 
assessments  of  AEO  from  measured  or  assessed  data  is  of  immense 
importance.  Worldwide,  wind  industry  is  trying  to  cope  with 
variance  in  AEO  assessments  from  the  actual  generation  achieved. 
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This  variance  can  be  reduced,  to  a  significant  extent,  by  paying 
greater  attention  to  the  frequency  distribution  models  assumed  in 
the  analysis. 

Therefore,  we  recommend  that  this  exercise  should  be  carried 
out  for  wind  speed  data  in  other  parts  of  the  world  and  additional 
formulations  established  for  undertaking  different  kinds  of  analy¬ 
sis  with  wind  speeds,  height  extrapolation,  wind  turbine  design 
considerations  and  energy  assessments. 
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